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Abstract 

The electronic structure of Co/Cu(001) and Fe/Cr(001) magnetic multi- 
layers has been investigated within the local density approximation com- 
bined with dynamical mean field theory. Our calculation shows enhanced 
density of states at the Fermi level, suggesting that electronic correlations 
might play an important role in the transport properties of multilayers. 
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I. INTRODUCTION 



Magnetic multilayers (MML) heterostructures of alternating ferromagnetic layers 
and non-magnetic spacers have attracted attention in the last decades because of their 
implications both for fundamental research and technological applications. The most 
remarkable property of these systems is the giant magnetoresistance (GMR) measured 
for a parallel/antiparallel configuration of the magnetic moments belonging to different 
layers in the presence of an external magnetic filed [1]. The GMR effect is related to 
the spin dependent scattering but the detailed mechanism is still subject of intense 
investigations. 

One of the first attempts to explain the GMR effect was based on the ballistic 
approach [2]. Soon it became obvious that the band structure description is very im- 
portant for the realistic description of GMR. A considerable number of attempts have 
been made to include the electronic structure in diffusive regime [3,4]. In order to con- 
sider the spin dependent scattering due to the interface, these calculations were based 
on the coherent potential approximation. However, all these calculations overestimate 
the experimental values for the GMR effect. On the other hand band structure cal- 
culations allows the evaluation of the Fermi velocities which can be used to estimate 
the spin dependent relaxation-times [5]. The semiclassical approach used in the cal- 
culation of the conductivity overestimate also, the experimental GMR values [6]. One 
of the most advanced technique in the ab initio theories of electric transport in solid 
systems, is offered by the Kubo-Greenwood formalism. A comprehensive overview of 
this technique applied for systems with reduced dimension is discussed by Weinberger 
[7]- 

In this paper, we demonstrate that the correlations effects might play an important 
role in the realistic description of transport properties. Our approach is motivated by 
the fact that the 3d transition metal elements, components of the magnetic multilayers, 
show significant correlation effects. In general, spin-polarized LDA band structure 
calculation gives an adequate description of the ferromagnetic ground state for the most 
of metals but, on the other hand, there are obvious evidences of essentially many-body 
features in photoemission spectra of Fe [8], Co [9], and Ni [10], [11]. The 6 eV satellite in 
Ni, broadening of the ARPES features due to quasiparticle damping, narrowing of the 
c?-band, essential change of spin polarization near the Fermi level are some examples. 
It is obvious that these effects can be equally important also for magnetic multilayers 
and other heterostructures containing transition metals. 

II. NUMERICAL DETAILS 

We focus on the band structure of the perfect (001) fee Coq/Cu^/Co^ and (001) bec 
Fe3/Cr 5 supercells. It was found experimentally that values of GMR are 220% in Fe/Cr 
multilayers [12] and 120% in Co/Cu multilayers [13]. One of the reasons why the above 
multilayers are highly magnetorezistve is that they contain ferromagnetic 3d metals 
which should have a pronounced spin asymmetry in their conductivity due to the pres- 
ence of the exchange split rf-bands. Perhaps the crucial factors for obtaining high values 
of GMR are the band matching and the lattice matching between the ferromagnetic 
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and non- magnetic metals [14,6]. These two conditions are almost perfectly satisfied 
in Co/Cu and Fe/Cr multilayers. Thin films of Co grow in the fee structure with the 
lattice parameter of 3.56A, which is only 2% less than the lattice parameter of 3.61% in 
fee Cu. Both Fe and Cr have the bec structure and their lattice parameters are almost 
identical: 2.87A for Fe and 2.77A for Cr. On the other hand it was experimentally 
demonstrated that the magnetic multilayers grown epitaxially shows an enhancement 
of the electronic contribution to the low temperature specific heat [15] which can be an 
evidence of correlation effects. Indeed, this enhancement cannot be reproduced in the 
standard electronic band structure calculations [16], probably because of the neglecting 
of many-body effects. It is our aim to check whether the correlation effects considered 
in the Dynamical Mean Field Theory (DMFT) [17] approach can lead to an essential 
renormalization of the density of state at the Fermi level N(Ep). 

We performed calculations based on the LDA+DMFT scheme [18]. The present 
realistic LDA+DMFT method [19] is based on the so-called Exact Muffin-tin orbitals 
(EMTO) method [20,21] within a screened KKR [22,23], frozen core together with 
the selfconsistent local spin density approximation (LDA). The correlation effects are 
treated in the framework of DMFT [17], with the Spin- Polarized T- matrix plus the Fluc- 
tuation Exchange (T-FLEX) approximation for the quantum impurity solver [24,25]. 

In the calculations for the Co/ Cu system we considered a tetragonal supercell formed 
by 16 layers, each layer containg one atom, with an interlayer distance corresponding 
to the fee Cu lattice constant 3.61A The supercell structure used for the study of the 
Fes/Cr^ system is presented in Fig. 1. As can be seen the optimization (relaxation) of 
the atomic layers are neglected, the interlayer distances correspond to the value of the 
bulk bec Fe 2.88A Each atom type is located on one layer, magnetically symmetric 
atoms are represented by the same colored spheres. Fe\ atoms, or interface atoms, are 
denoted by blue sphere. The central, Fe 2 atoms are indicated by a light blue color. 
In the picture representing the structure Fig. 1 three extra Fe layers: Fe\j Fe^j Fe\ 
belonging to the next unit cell were introduced. Interface CV li2 layers are indicated by 
green, respectively yellow color spheres, meanwhile the central CV 3 layer is denoted by 
red sphere. In the calculations the same atomic sphere radius was considered for Fe 
and Cr in the Fe/Cr, and similarly the radii for Co and Cu atoms were chosen the same 
for the Co/Cu system. 

The 3s, 3p, and 3d states are included for the Co, Cu and Fe valence electrons. In 
order to calculate the charge density we integrate the Green function along a contour on 
the complex energy plane which extends from the bottom of the band up to the Fermi 
level [21], using 30 energy points. For the Brillouin zone integration we sum up a number 
of 567 k-points for the Co/Cu system and 839 k points for the Fe/Cr respectively. A 
cutoff of Imax — 8 for the multipole expansion of the charge density and a cutoff of l max = 
4 for the wave function was used. The Perdew-Wand [26] parameterization of the local 
density approximation to the exchange correlation potential was used. Although the 
lattice mismatch and relaxation have been shown to influence the magnetic properties 
[27] we neglect this effect. 
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A. fee Co and Co/Cu multilayers 



We will first discuss the correlation effects on bulk fee Co, and afterwards investigate 
the Co/Cu multilayer. The LSD A and LSDA+DMFT density of states for fee Co is 
shown in Fig. 2. The LSDA electronic structure of fee cobalt is different for the majority 
and minority spin electrons, due to it's ferromagnetism \ilsda = 1.69/ie. DOS can be 
characterized qualitatively by a shift of the minority and majority c?-bands relative 
to each other. Due to this asymmetry the contributions of the two spin channels to 
the density of states at the Fermi level will be different, therefore the conductivities 
are different. This asymmetry in the DOS is the source of GMR in magnetic Co/Cu 
multilayers. The most important feature of the cobalt DOS is that the Fermi level lies 
above the top of the ci-band for the majority spin electrons. 

For the LSDA+DMFT calculations a value of the average Coulomb interaction 
U — 2 eV was chosen. Correlation effects related with spin-flip excitations at non- 
zero temperature reduce the magnetic moment, ^dmft = l-42/z B in comparison with 
the corresponding LSDA value. The asymmetry between the majority /minority spin 
channels is kept and an increasing of 12% of DOS at the Fermi level N(Ep) is evidenced. 
At low energies the majority spin channel presents a satellite structure similar to the 
one found in the Ni. Due to the perturbation nature of our approach the satellite is 
shifted towards lower energies. 

The energy dependence of the self-energy for Co is plotted in Fig. 3. Near the Fermi 
level a typical Fermi liquid behavior is evidenced. For the imaginary part —Im S(-E') oc 
E 2 , meanwhile the real part of the self-energy has a negative slope dRe H(E)/dE < 0, 
where E is the electron energy relative to the Fermi level. Due to the fee structure 
Co self-energy shows a considerable similarity to the selfenergies of Ni [19]. The high 
value of the imaginary part of self-energy evidenced around — 7eV in the majority spin 
channel for both t2 9 and e g orbitals produce the satellite visible in DOS Fig. 2. 

The Coq/ Cu§/ Co?, multilayer, total and the interface layers DOS, in the ferromag- 
netic orientation are presented in Fig. 4. As we seen in Fig. 4, the DOS is asymmetric 
between the majority and minority spins. In comparison with the bulk fee Co , Fig. 2, 
a qualitative description of DOS of Co 6 /Cu 5 /Co 5 can be made, on the base of a "rigid" 
shift due to the presence of the Cu layers. The Co layers magnetic moment increase 
as we approach the interface, meanwhile the Cu layers are non-magnetic; the values of 
the magnetic moments and the electronic specific heat coeficients are presented in table 
I. The latter is given by the relation 7 = ir 2 k 2 B N(E F )(l + A)/3 where N(E F ) is the 
electronic DOS at the Fermi level, (1 + A) is the mass enhancement factor caused by the 
electron-phonon interaction. The calculated LSDA+DMFT magnetic moments in the 
Coq/Cus/Cos superlattice is shown in Fig. 6, where for comparison we also plot the 
corresponding quantities obtained using the LSDA. We note that at finite temperatures 
the Cu layers remain non-magnetic, whereas the value of the magnetic moment on the 
Co layers diminish in comparison with the LSDA one. Even at finite temperatures 
(250^) the Co moments couples ferromagnetically across the Co/Cu interfaces. 

The energy dependence of the self-energy for the Co interface and central layers are 
plotted in Fig. 5. Near the Fermi level a typical Fermi liquid behavior is evidenced 
for both the interface and central Co layers. The imaginary part — Im T,(E) oc E 2 , 
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meanwhile the real part of the self-energy has a negative slope dRe / BE < 0, where 
E is the electron energy relative to the Fermi level. As can be seen the selfenergies and 
implicitly the correlation effects are different for the distinct Co layers. 

B. bcc Fe, bcc Cr and Fe/Cr multilayers 

The correlation effects on bulk bcc Fe, bcc Cr and (OOf ) surface of bcc Fe covered 
by a trilayer of Cr were described in a previous paper [19]. The present calculations 
addresses the finite temperature and correlations effects in an eight-layer (3Fe+5Cr) 
bcc superlattice. 

In Fig. 7 the bcc bulk Fe and Cr DOS are presented. As we can see due to the 
bcc structure the DOS exhibits a pronounced valley in the middle of the ci-bands for 
both spins. Fermi level lies within the ci-band for both spin orientations which provides 
dominance of the (^-character. Iron is magnetic Hlda = 2.25/xb, whereas chromium is 
non- magnetic. 

Similarly to the bulk DOS, as it is evident from Fig. 8 the Fermi level in the Fe/Cr 
system lies within the ci-band for both spin orientations. DOS exhibits a pronounced 
valley for the minority spins with the Fermi level lying almost at the bottom of this 
valley. This feature of the band structure is a consequence of the similarity of the DOS 
of minority spin electrons of iron and the minority spin channel DOS of chromium Fig. 
7. 

Many-body effects described in the framework of the LDA+DMFT, were investi- 
gated for different temperatures: T = 250, 500i^. Since we expect similar correlation 
effects in the Fe/Cr system the same value of the average Coulomb interaction U = 2eV 
and the same exchange correlation energy J = 0.9eV was chosen. The DOS for the Fe 
layers and Cr layers are presented in Fig. 8 and Fig. 10 respectively. 

In the case of Fe, Fig. 8 the correlation effects are manifested in a slightly different 
way for central and interface Fe layers. One can see that for the interface (Fei) layer 
the LSDA peack of the unoccupied DOS close to the Fermi level, is pinned at the Fermi 
level producing an enhancement of N(Ep). The spectral weight of the main peak in 
the occupied part (— leV) of spin up channel, is transfered closer to the Fermi level, as 
temperature is increased. In the same time finite temperature effects smear out the low 
energy features of DOS, situated in the energy range of —4, —2 eV. In the spin down 
channel, as the correlation effects are switched on a peak at the Fermi level appears. 
The spectral weight of the leV spin down LSDA peak is transfered towards the Fermi 
level and a new peak appears around 0.25eV. At T = 250K, Fig. 8, the weight of 
the former leV peak dominates the weight of the 0.25eV peak, but as the temperature 
increases, T = 500i^, more spectral weight is transfered to the 0.25eV peak. 

Similarly to the interface Fe\ layer, the spin down density of states of the central 
Fe2 layer show the appearance of the peak at the Fermi level, and the temperature 
dependence of the spectral weight transfer towards the Fermi level. Due to the corre- 
lation effects, the narrowing of the width of the —leV peak in the spin up channel is 
evidenced. As the temperature is increased, the — leV peak shows a slight shift towards 
the Fermi energy, but it's width is not significantly changed. 
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The energy dependence of the self-energy for the Fe layers are plotted in Fig. 9. 
Near the Fermi level a typical Fermi liquid behavior is evidenced for both the interface 
and central Fe layers. The imaginary part —Im S(-E') oc E 2 , meanwhile the real part of 
the self-energy has a negative slope dRe H{E)/dE < 0, where E is the electron energy 
relative to the Fermi level. As can be seen the selfenergies and implicitly the correlation 
effects are different for the two distinct Fe layers. 

All the above correlation effects can be recognized in the Cr layers as well Fig. 10, 
in particular the formation of the peak at the Fermi level being a significant feature 
for the spin down channel. According to the LSDA calculation the density of states 
of the central Cr layer (Cr 3 ) is almost the same as that of the bcc Cr bulk Fig. 7. As 
going from the interface Cr layers (CVi and Cr-i) towards the central layer (CV3) the 
spin up channel density of states show the formation of a valley in the DOS near the 
Fermi level characteristic to the bulk behavior Fig. 10. The interface, Cr\ layer, DOS 
is modified appreciable because of the presence of nearby strongly ferromagnetic Fe 
layer. It is important to note, however that the spin down channel of all Cr layers 
are not significantly affected by the proximity to the magnetic Fe layer. The spectral 
weight transfer is also present, but for the interface Cr\ layer is not that evident since 
hybridization with the nearby strongly ferromagnetic Fe layer is present. 

In Fig. 12 we notice that the magnetic moments of Cr layers alternate from layer 
to layer and the Fe moments couple ferromagnetically across the Fe/Cr interfaces. 
According to our LDA+DMFT calculation this Fe/Cr antiparallel coupling across the 
interface was estimated to be more stable compared to the parallel solution. The 
temperature dependence of the Fe magnetic moment on each layer is presented in 
Fig. 13. The central Fe<i layer follows approximately the behavior of the bulk values, 
meanwhile the interface layer Fe\ has a faster temperature decrease. On the other 
hand the Cr layers show a very peculiar temperature dependence of the magnetic 
moment due to the couplings across the Fe/Cr layers. Fig. 14 display the Cr layers 
temperature dependence of the magnetic moment which shows a deviation from the 
Brillouin function. 

The Cr layers being non-magnetic the majority and minority spin self-energies 
Fig. 11 are identical. However correlation effects seems to be more important for the 
central CV 3 layer. 

III. CONCLUSIONS 

The mechanism of giant magnetoresistence in magnetic multilayers is usually re- 
lated to the spin- dependence of the scattering process. The spin dependent scattering 
is assumed to arise from spin-dependent random potentials produced by magnetic im- 
purities at the interface or in the bulk of the ferromagnetic layers [34,35]. Recently 
an improved prediction of the GMR was obtained by combining the disorder effects 
with the accurate spin-dependent electronic structure calculations [14]. However finite 
temperature properties and correlations effects were not taken into account. 

Several theoretical approaches have been used to explain the magnetic properties 
of such superlattice structures. Many of these approaches are based on the RKKY-like 
model [31,32], tight-binding models [33] and, recently, on the results of ab initio elec- 
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tronic structure calculations [29]. The magnetic coupling studied in the framework of 
these models was shown to result from the interplay between the direct d — d hybridiza- 
tion of Fe and Cr atoms and indirect exchange through the sp electrons. The sp — d 
coupling [29] was found to be reminiscent of the RKKY interaction only for superlat- 
tices with more than four Cr layers [29]. However, for the case of three Cr layers the 
sp — d coupling model cannot explain the ferromagnetic ordering of Fe atoms. 

Based on the theory of spin-fluctuations developed by Moriya [38] and applied by 
Hubbard [39] for the case of several transition metals, Hasegawa [37] showed that spin- 
fluctuations plays an important role in discussing the temperature dependence of the 
GMR. Due to the static approximation employed in the model [37], dynamical effects 
of the spin-fluctuations were neglected, therefore any definite conclusion on the tem- 
perature dependence of GMR were not drawn. Although the explicit calculation of 
the resistivity and GMR is not the purpose of the present paper, our results include 
dynamical correlations described in the framework of DMFT [17], being a promising 
starting point for a direct evaluation of the GMR in multilayer systems. 

Using an first principle LDA+DMFT approach [19] we examined the correlation 
effects and the finite temperature magnetic properties of some Co/Cu and Fe/Cr 
superlatticels. Our calculations evidenced a peculiar temperature dependence of the 
magnetic moments near the interface. The correlation effects proved to be different for 
different atomic layers. Based on the spin polarized T-matrix FLEX, DMFT solver, 
the correlation effects capture the spin fluctuations that plays primary roles at finite 
temperatures. Recent experiments on Fe/Cr trilayers showed the existence of magnetic 
fluctuations [36]. Even though the magnetic excitation were attributed to structural 
and magnetic disorder in the vicinity of both Fe/Cr interface [36], results presented in 
this paper suggest that electron-electron interactions give rise to magnetic excitations, 
which are common features in 3d transition metals systems. 

It is worthwhile to emphasize that the enhanced DOS at the Fermi level, having a 
many-body correlation origin, can play an important role in the GMR, since this DOS 
enhancement is strongly spin-dependent. It is more effective in the minority channels 
of Co/Cu And Fe/Cr systems, giving the result of a quasiparticle peak centered at the 
Fermi level. Our calculation is in good agreement with the tendency of the enhancement 
of electronic contribution in Fe/Cr magnetic multilayers [15]. 
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TABLES 



TABLE I. Theoretical magnetic moments fx and electronic specific heat coeflcients 7 cal- 
culated for bulk bcc Fe, Cr and fee Co. The experimental values for 7 include the contribution 
from the electron-phonon coupling, which is not included in the results obtained from band 
structure calculations. In the last three columns, the parameters used in the self-consistent 
LSDA+DMFT calculations are listed. 





I^LDA 


UDMFT 


1LDA 


1DMFT 


7errp 


T 


U 


J 




(Mb) 


(vb) 


(mJ/K 2 mol) 


(mJ/K 2 mol) 


(mJ/K 2 mol) 


(K) 


(eV) 


(eV) 


Co 


1.62 


1.38 


5.43 


6.78 




250 


2 


0.9 


Fe 


2.25 


2.24 


2.43 


2.61 


3.11,3.69 a 


250 


2 


0.9 


Cr 






1.73 


2.40 


3.5 b 


250 


2 


0.9 



a Ref. [40] 

b Non- magnetic Cr, Ref. [41]. 



TABLE II. 


The calculated 


magnetic moments in the LSDA (T=0K) and LSDA+DMFT 


(T=250K) for 


the Coq/Cu^/C 


05 superlatice. 








I^LDA 


Vdmft(T = 250K) 


U 


J 




(vb) 


(p-b) 


(eV) 


(eV) 


C 0l (Co/Cu) 


1.44 


1.12 


2.0 


0.9 


Co 2 (Co/Cu) 


1.44 


1.02 


2.0 


0.9 


Co 3 (Co/Cu) 


1.37 


1.04 


2.0 


0.9 


Co A {Co/Cu) 


1.31 


1.06 


2.0 


0.9 


Co 5 (Co/Cu) 


1.31 


1.09 


2.0 


0.9 


Co 6 (Co/Cu) 


1.35 


1.20 


2.0 


0.9 


Cui{Co/Cu) 


0.00 


0.02 


0.0 


0.0 


Cu 2 {Co/Cu) 


0.00 


0.00 


0.0 


0.0 


Cu 3 (Co/Cu) 


0.00 


0.00 


0.0 


0.0 



TABLE III. Temperature dependence of layer resolved magnetic moments. The bulk 
values of the magnetic moment for the same temperatures, and the parameters used in DMFT 
calculations are presented in the table. 





VLDA 

(Mb) 


li{T = 250K) 
(Mb) 


M (T = 500K) 
(Mb) 


U 
(eV) 


J 

(eV) 


Fe(bulk) 


2.25 


2.24 


2.21 


2.0 


0.9 


Fei{Fe/Cr) 


1.97 


1.78 


1.70 


2.0 


0.9 


Fe 2 (Fe/Cr) 


2.54 


2.38 


2.42 


2.0 


0.9 


Cri(Fe/Cr) 


-0.54 


-0.04 


-0.08 


2.0 


0.9 


Cr 2 (Fe/Cr) 


0.43 


0.05 


0.04 


2.0 


0.9 


Cr 3 (Fe/Cr) 


-0.45 


-0.18 


-0.16 


2.0 


0.9 
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FIGURES 




FIG. 1. The superlattice strcuture of Fe^/Cr^ having an equidistant distribution of atomic 
layers. Fe atoms are denoted by blue spheres and Co atoms are indicated by green yellow 
and red sphere, in the following order: Fe\jFe<il 'Fei/Cri/Cr2/Crs/Cr2/Cri/ 'Fe\/ 'Fe^jFex. 
The tetragonal supercell is formed by the first 8 atomic layers (3Fe+5Cr) aligned along the z 
direction. 




FIG. 2. The LSDA (dashed line) and LSDA+DMFT (solid line) densities of states for 
fee Co calculated using the EMTO-DMFT method. A significant reduction of the exchange 
splitting can be evidenced. 
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FIG. 5. Layer resolved self-energies: C 06-interface layer and Co\ - central layer at tem- 
perature T = 250. 
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FIG. 6. Distribution of the layer resolved magnetic moment at T=0K (LSDA) and 
T=250K (LSDA+DMFT) in Co/Cu super lattice. 
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FIG. 8. Total DOS for both types of Fe layers: Fei-interface layer and Fe2 - central layer. 
The temperature dependence of DOS is presented for T = 0, 250 and 500K respectively. 
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FIG. 9. Layer resolved self-energies: Fei-interface layer and Fe2- central layer at temper- 
ature T = 500. 
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FIG. 10. Total DOS for the three types of Cr layers: Cr±,Cr2 - interface layers and CV3 
- central layer. The temperature dependence of DOS is presented for T = 0, 250 and 500K 
respectively. 
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FIG. 11. Layer resolved self-energies for Cr layers: Cr\,Cr2 - interface layers and CV3 
central layer for for T = 500 . 
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FIG. 12. The distribution of the local magnetic moments at T=0, 250 and 500K in 3Fe 
+5 Cr superlattice as given by the LSDA and the different DMFT calculations 
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FIG. 13. Temperature dependence of the calculated magnetic moments on the Fe layers is 
presented in comparison the temperature dependence of the calcualted bulk bcc Fe magnetic 
moments. The value at T = 325-ftT is taken from the reference [19]. The normalization of the 
magnetic moment is done with respect to the LSDA value (/j>t=o)- 
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FIG. 14. Temperature dependence of the magnetic moment on the Cr layers. The nor- 
malization of the magnetic moment is done with respect to the LSDA value (jPr=o). 
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